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ABSTRACT 

We present a detailed structural and morphological study of a large sample of spectroscopically- 
confirmed galaxies at z > 6, using deep HST near-IR broad-band images and Subaru optical narrow- 
band images. The galaxy sample consists of 51 Lya emitters (LAEs) at z ~ 5.7, 6.5, and 7.0, and 16 
Lyman-break galaxies (LBGs) at 5.9 < z < 6.5. These galaxies exhibit a wide range of rest-frame UV 
continuum morphology in the HST images, from compact features to multiple component systems. 
The fraction of merging/interacting galaxies reaches 40 — 50% at M1500 ^ —20.5 mag, suggesting 
hierarchical build-up of the brightest galaxies at high redshift. We use half-light radius to describe the 
galaxy sizes, and find that the intrinsic radii rhi^n, after correction for PSF broadening, are roughly 
between rhi,in — 0'.'05 and O'.'S, with a median value of O'.'IS (~0.75 kpc). This is consistent with the 
sizes of bright LAEs and LBGs at z > 6 in previous studies. In addition, more luminous galaxies 
tend to have larger physical sizes, exhibiting a size-luminosity relation rhi^n oc L°'^. The slope 0.2 
is significantly fiatter than those in previous fainter LBG samples. We for the first time characterize 
the morphology of z > 6 galaxies using nonparametric methods, including the CAS system, the Gini 
and M20 parameters. Compared to low-redshift galaxies, our galaxies appear in slightly different 
locations in the parameter space, mainly due to their small sizes in the HST images. However, we 
find strong correlations between the measured parameters, as expected from their definitions. This 
implies that these nonparametric methods could be still applicable for z > 6 galaxies, if used with 
caution. We search for extended Lya emission halos around LAEs at z ~ 5.7 and 6.5, by stacking a 
number of narrow-band images. We do not find evidence of extended halos predicted by cosmological 
simulations. Such Lya halos, if exist, could be weaker than predicted. Finally, we investigate any 
positional misalignment between UV continuum and Lya emission in LAEs. While the two positions 
are generally consistent, several merging galaxies show significant positional differences. It is likely 
caused by disturbed ISM distribution due to merging activity. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: high-redshift 
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Galaxy structural and morphological studies provide 
basic apparent information about galaxies. Nearby 
galaxies are generally classified into three broad cate- 
gories, spiral, elliptical, and irregular galaxies. The ma- 
jority of lumino us nearby galaxies (z < 0.1 ) are spirals 
and ellipticals ([Abraham fc van den Ber"ghl |2001[ ). At 
higher redshift, galaxies are less well developed, and 
the fraction of irre gular galaxies increases steadily (e.g., 
[Driver et ai]|1995l 1X998'). In the redshift range of 0.1 < 
z < 1, galaxy morphology and structu re have been well 
studied (e.g. jBrinchmann ct al. 199 ^ iLillv et al. 199| 
Schade et aL|[l99a: iCarlberg et al... 20001: iLe Fevre et al 



2000n van den Bergh et al.ll2000l . and references therein) 



These galaxies show more disturbed structures than 
nearby galaxies do in the rest-frame UV and optical (e.g 



Abraham fc van den Bcrgh'200r: Windho rst et al, 

Blanton fc Moustakai 



2002 



2008; 



Tavlor-Mager et all ,2007, 

Shi et al.l |2009() . The fraction of irregular galaxies in- 
creases from less than 10% at z < 0.5 to ~ 30% at 



z^l (e.g.. lBrinchmann et al.|[T998l : van den Bergh et al.l 
l2003). In addition, more galaxies were identified as merg- 
ing systems, refiecting hierarchical build-up of galaxies 
and mass assembly in cold dark matter (CDM) scenario 
(|White fc ReesI [19781 iCole et al.l [2000l). For example 
nearly 20% of z ~ 1 galaxies in the ILe Fevre et al.l ([2000l ) 
sample are in close pairs. 
For galaxies at z > 2 — 3, morphological classifica- 
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tion is challenging, as most galaxies appear peculiar. 
Galajcies are also s maller towards higher redshift 
(|Ferguson et al.l I2004D . In addition, galaxies appear 
much fainter due to the cosmological (1 + z)'^ sur- 
face brightness dimming. Traditional classifications, 
including Hubble's tuning-fork system, do not work 
any more. Theref ore, nonparamet ric methods such as 
the CAS system (iConselice 120031) and the Gini and 
M20 parameters (jLotz et al.l 120041 ) play an important 
role. Most morphological and structural analyses 
in this redshift range were done in the GOODS fields 
(jGiavalisco et al.ir2004|) , because of the high-quality Hub- 
ble S pace Telescop e (HST) d ata (e .g .. iLowenthal et aJ] 
iRavindran ath et all 120061 : iLaw et al.l 120071 : 



T997I: 



Cassata et al.l l201^ In the rest-frame UV, z > 3 



galaxies are usually compact (from one to several kpc), 
but many of them display extended fea tures or multiple 
clumps in deep HST images (e. g., [G iavalisc o et al. | 
1996 ; iVenemans et all l2005t IRavin dranat h et al.l 
2006t IPh^al et al.l I2OO7I IConselice fc Arnold l200a 
Gronwall et al, iLaw et al ji2oia). For example, the 

Ravindranath et al.l ( 2006D sample contains thousands 



of photometrically-selected (hereafter photo-selected) 
LBGs at 2; > 2.5, and abou t 30% of them have 
multiple cores. In I Law et aTl (|2012l ). a sample of 
spectroscopically-confirmed (hereafter spec-confirmed) 
galaxies also showed a high fraction of interacting 
systems. 

In the highest-redshift range z > 6, morphological 
studies become difficult. Galaxies appear very faint, and 
their rest-frame UV light moves to the near-IR wave- 
length range, where telescope resolution is poorer. A 
typical galaxy occupies only a few pixels even in HST 
near-IR images, so nonparametric parameters have not 
been determined for these galaxies, and size is usually 
the only physical parameter that can be reliably mea- 
sured in the literature. Studies based on photo-selected 
galaxies have shown that z > 6 galaxies are generally 
very compact, and most of them arc just barely spatially 
resolved. For example, Ocsch et al. (2010) reported the 
sizes of 16 LBGs at z > 7 in the Hubble Ultra-Deep 
Field (HUDF). They found that only two in their sample 
show extended features, and the rest are very compact 
(< 1 kpc). This sample is very faint. Observations of 
a handful brighter galaxies with spectroscopic redshifts 
also suggest c ompact morphology, with a typical size of 
< 1 kpc fe.g..lStanwav et al.ir2CI04l: iDow-Hvgelund et al.l 
120071: ICowie et al.l 120 lit) . Note that galaxy size IS corre- 
lated with physical properties, such as mass and lumi- 
nosity at low redshift. Such relations still exist in high- 
redshift galaxies, but could have evolved with time (e.g., 
iGrazian ct al. 2012; Moslch et al. 2012). 

In this paper, we will carry out a structural and mor- 
phological study of a sample of 67 spec-confirmed galax- 
ies at z > 6. The sample is the largest collection of 
spec-confirmed galaxies in this redshift range, including 
51 Lya emitters (LAEs) and 16 Lyman-break galaxies 
(LBGs). This paper is the second in a series present- 
ing the physical prop erties of these gal axies. In the first 
paper of the series (jJiang et al.l 120131 hereafter Paper 
I), we presented deep Subaru optical and HST near-IR 
data. We also derived various rest-frame UV continuum 
and Lya emission properties, including UV-continuum 



slope /3, the Lya rest- frame equivalent width (EW), and 
star formation rates (SFRs). These galaxies have steep 
UV continuum slopes roughly between /? ~ —1.5 and 
(3 ~ —3.5, with a median value of /? ~ —2.3. They have 
a range of Lya EW from --10 to --300 A. Their SFRs are 
moderate from a few to a few tens solar masses per year. 
In this paper we will study the structure and morphology 
of their rest-frame UV-continuum emission based on our 
HST images, and of their Lya emission based on our 
ground-based narrow-band images. 

The layout of the paper is as follows. In Section 2, 
we briefiy review our galaxy sample and the optical and 
near-IR data that will be used for the paper. We measure 
the structure and morphology of UV continuum emission 
in Section 3, and of the Lya emission in Section 4. We 
then discuss our results and summarize the paper in Sec- 
tion 5. Throughout the paper we use a A-dominated fiat 
cosmology with Hn = 70 km s~ ^ Mpc~^, ftm = 0.3, and 
r^A = 0.7 (iKomatsu et al.l 120111). A ll magnitudes are on 
the AB svstem (|Oke fc Gunnlll9"83l ). 



2. GALAXY SAMPLE AND DATA 

In Section 2 of Paper I we provided a detailed descrip- 
tion of our galaxy sample, and the data that we used in 
multi-wavelength bands. Here we summarize this infor- 
mation below. There are a total of 67 spec-confirmed 
galaxies in our sa mple: 62 are fro m the Subaru Deep 
Field (SDF; Kashi kawa et al.ll2004[ ). and the remammg 
5 are from the Subaru XMM-Newton Deep Survey field 
(SXDS: .Furusawa et al.,, 2008 1. They represent the most 
luminous galaxies in terms of Lya luminosity (for LAEs) 
or UV continuum luminosity (for LBGs) in this redshift 
rang e. The SDF galaxy sample contains 22 LAEs at z ~ 
5.7 (jShimasaku ct al. 2006; Kashikawa et al. 201 1), 25 
LAEs at z ~ 6.5 (Taniguchi ct al. 2005; Kashikawa ct gJ] 
I2OO6I [201lh . and a LAE at z = 6.96 (|lve et al.l [2001 7 
The LAEs at z ~ 5.7 and 6.5 have a relatively uni- 
form magnitude limit of 26 mag in the narrow bands 
NB816 and NB921, and thus make a well-defined sam- 
ple. Th e SDF sample also conta ins 14 LBGs at 5.9 < 
z < 6.5 (INagao et a l.''2004. '2005, 2007: Ota et al.l 120081 : 
iJiang et al.l 120111 : iToshikawa et al . 20121 These LBGs 
were selected with different criteria, and have a rather 
inhomogeneous depth. The SXDS sample consis ts of 5 
galax ies, including 2 LBGs at z :^ 6 (Curt is- Lak e et al.l 
120 12f ) and 3 LAEs at z ~ 6.5 (jOuchi et al.ll201Q) . In Ta- 
ble 1 of Paper I, we summarized all the observed data of 
these galaxies. 

The SDF a nd SXDS were o bserv e d with Subaru 
Suprime-Cam (jKashikawa et al.l 120041 : iFurusawa et al.l 
They have extremely deep optical images in a 



series of broad and narrow bands. Public stacked images 
are also available for the two fields, but the public data 
do not include the images taken recently. In Paper I, 
we produced our own stacked images in six broad bands 
(BVRi'z'y) and three narrow bands (NB816, NB921, 
and NB973) by including all available data in the archive. 
Our stacked images have great depth with excellent PSF 
full width at half maximum (FWHM) of 0'.'5 - 0'.'7. The 
near-IR imaging data for the SDF galaxies are from three 
HST GO programs 11149 (PI: E. Egami), 12329 and 
12616 (PI: L. Jiang). The HST observations were made 
with a mix of instruments and depth. The majority of the 
galaxies were observed with WFC3 in the F125W (here- 



3 



after J125) and F160W (hereafter H or -ffieo) bands. The 
typical integration time was two HST orbits (roughly 
5400 sec) per band. This provides a depth of ~ 27.5 mag 
(5a detection) in th e J12R band and ^ 27.1 mag in the 
HiQQ band (see also [Windhorst et al.l 120 111) . The pixel 
size in the final reduced WFC3 images is 0'.'06. Sev- 
eral SDF galaxies were observed with NICMOS in the 
FllOW (hereafter Jno) and iJieo bands. The typical in- 
tegration time was also two HST orbits. The depth in 
the two bands are ~ 26.4 mag and ~ 26.1 mag, respec- 
tively. The pixel size in the final reduced NICMOS im- 
ages is O'.'l. The five SXDS galaxies were covered by the 
UKIDSS Ultra-Deep Survey (UDS). Their HST WFC3 
near-IR data were obtained from the Cosmic Assembly 
Near-i nfrared Deep Extrae alactic Legacy Survey (CAN- 
DELS; iGrogin et aT1l201lh . The exposure depth of the 
CANDELS UDS data is 1900 sec in the J125 band and 
3300 sec in the Hiqq band. 

The majority of the galaxies in our sample were de- 
tected with high significance in the near-IR images. Only 
15 of them — among the faintest in the optical — have 
weak detections (< 5a) in the J band (J125 or Jno)- 
In Table 1 of Paper I we listed the optical and near-IR 
photometry of the galaxies. In Table 2 of Paper I we pre- 
sented basic physical properties, including the rest-frame 
UV continuum luminosity and slope /3, the Lya luminos- 
ity and EW, and their SFRs, etc. The thumbnail images 
of all the galaxies are provided in Appendix A of Paper 
1. 

3. UV CONTINUUM MORPHOLOGY 

In this section we will derive structural and morpho- 
logical parameters for the galaxies in our sample. Al- 
though our galaxies represent the most luminous galax- 
ies at 2 > 6, they appear faint and small compared to 
lower-redshift galaxies. The majority of them are point- 
like sources in the optical images. Even in the HST 
WFC3 images they usually occupy a very limited num- 
ber of pixels. Hence, the study of these distant galaxies is 
challenging. In previous literature, galaxy size was often 
the only parameter that could be reliably measured for 
z > 6 galaxies. Classifications for nearby galaxies, such 
as the classical Hubble's tuning-fork system, cannot be 
applied to these objects. In this section, we will measure 
the sizes of our galaxies, and try to characterize their 
morphology using nonparametric methods, such as the 
CAS syst em LCo nselice 2003) and the Gini and M20 pa- 
rameters ()Lotz et al.ll2004f ) . These methods are primarily 
used for low-redshift galaxies, though they have already 
been used for galaxies at 2: = 4 ~ 6 (e.g., IPirzkal et al.l 
I2007t IConsehce fc Arnoldll2009H . We will also study in- 
teracting/merging systems in our sample. 

In order to calculate the above parameters, we take 
advantage of all our HST images. For each galaxy 
in our sample, we combine (the weighted average) its 
J- and i^-band images and make a stacked HST im- 
age to improve the signal-to- noise (S/N) ratio. By do- 
ing this we assume that the effect of the morphological 
k-correction — the dependence of galaxy structure on 
wavelen gth — is negligible in th e wavelength range con- 
sidered (j Windhorst et al.ll2002D . This is because the J 
and H bands cover a similar rest-frame UV wavelength 
range (~1780 A vs. ~2200 A) for 2; ~ 6 galaxies. Our 
further analyses are then based on the stacked images. 



The J and H bands do not cover Lya emission for our 
galaxies, so the morphology in the stacked images is 
purely from their UV continuum emission (other neb- 
ular lines can be safely ignored in general; see also Cai 
et al. 2011 and Kashikawa et al. 2012). In Figure 1, 
we show the thumbnail images of 47 (out of 67) galaxies 
that have more than 8ct detections of their total fluxes 
in the stacked images. We will focus on these 47 galaxies 
in this section. Table 1 lists the 47 galaxies. Columns 
1 to 4 list the galaxy number, coordinates, and redshift, 
directly taken from Paper I. Note that we exclude object 
no. 12, since it overlaps with a bright foreground star, 
as explained in Paper I. 

We do not show in Figure 1 the galaxies (19 out of 
67) that have < 8a detections, as we would not obtain 
reliable structural and morphological information with 
such low S/N ratios. Instead, we combine their images, 
and make a stacked image for these 19 faint galaxies. 
The individual images are scaled before stacking, so that 
the galaxies all have the same magnitude. Because of the 
scaling, the individual weight images have different sky 
noise. The final stacked image has a much higher S/N 
ratio. 

3.1. Size 

We use the half-light radius rhi to describe the size of 
a galaxy. The half-light radius rhi is a radius enclosing a 
half of the total light. The total light is measured within 
a radius of 1.5 rp, where rp is the Petrosian radius, at 
which the surface brightness (SB) is 20% of the aver- 
age SB within this radius. The factors 1.5 and 20% are 
chosen to ma tch those in previous studies of high-redshif t 
galaxies (e.g. lLotz et al.ll2006l: IConselice fc Arnoldll2009D . 
The Petrosian radii rp for the majority of the 47 galaxies 
in our sample are between 0'.'3 and 0'.'7, consistent with 
th e radius range 0"2 ^ 0'.'6 in the 2; ~ 4 — 6 LBG samples 
of IConselice fc Arnoldl (|2OO90 . 

Column 5 in Ta ble 1 lists the observed r hi calculated 
using SExtractor (jBertin fc ArnoutsI [19961) . The upper 
panel of Figure 2 shows ry as the function of M1500, the 
absolute AB magnitude of the continuum at rest-frame 
1500 A. Physical quantities such as A/1500 were derived in 
Paper I. The blue and red circles in Figure 2 represent the 
LAEs at z ~ 5.7 and 6.5 (including z ~ 7), respectively, 
and the green circles represent the LBGs at z ~ 6. These 
galaxies span a range of rhi from O'.'l to ^ 0'.'28, or 0.6- 
1.6 kpc, without correction for PSF broadening. The 
median value of the measured rhi is 0'.'17 1 kpc). The 
stacked image of faint galaxies, shown as a star in Figure 
2, has a relatively smaller radius rhi — 0'.'14. The median 
value of rhi for the whole sample is 0'.'16. 

We correct for PSF broadening and estimate intrin- 
sic half-light radius rhi, in in quadrature, i.e., r^ ;„ = 

'"hi ^ '"psfj where rpgp is the radius of the instrumen- 
tal PSF, and is O'.'IO in our HST WFC3 images. The 
above equation is a good approximation, if the PSF pro- 
file and the intrinsic object profile are Gaussian-like or 
other simple symmetric profiles. As we will see, some 
galaxies in our sample show rather asymmetric and ex- 
tended features. These galaxies usually have sizes much 
larger than the PSF, so their measured rhi are close to 
''hi in- With this correction, the intrinsic half-light radii 
ru,in range from < 0'.'05 (< 0.3 kpc) to - 0'.'26 (~ 1.5 
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kpc), with a median value of 0'.'14 (~ 0.8 kpc). The rhijn 
in the stacked image of 19 faint galaxies is about O'.'l. 
The median value of rhi.in for the whole sample is O'.'IS 
(~ 0.8 kpc). The results are plotted in the lower panel 
of Figure 2. 

The galaxy sizes in our sample roughly agree with those 
of high-redshift LAEs and LBGs with sim ilar luminosi- 
ties in the previous literature. For example, iPirzkal et al.l 
(|2007D found that the average r ^] for a sample o f lumi - 
nous LAEs at z ~ 5 is 0'.'17. iTaniguchi et al.l (|2009l ) 

fou nd a median rhi of O'.'IS for LAEs at z ^ 5 .7. In 

the iHathi et al.l pOOl and iConselice fc Arnoldl (poM ) 
LBG samples of z = 4 6 galaxies, the rhi range 
is o'.'l ^ 0'.'3, similar to the rhi range in our sam- 
ple. Previous studies have shown that the galaxy size 
roughly scale s with redshift a s (1 4- z)~"^, with m close to 
1.1-1.2 (e.g . iFerguson ct al.' 2004; B ouwens et~all 120061: 
lOesch et all [2010; Mosleh ct_al. 20 ll), so the size of 
galaxies evolves slowly at high redshift. This is the 
reason that high-redshift g alaxies have a similar size 
range. iMalhotra et al.l (|2012f ) found, however, that LAEs 
have a roughly constant size in the redshift range of 
2.25 < z < 6, and do not show a size- redshift relation. 
While our sample does not have a large redshift cover- 
age, our galaxy sizes are well consistent with those in 
their sample. 

3. LI. Size-luminosity relation 

Figure 2 shows that brighter objects tend to have larger 
Thi, meaning that more luminous galaxies have larger 
physical sizes. This size-luminosity relation has been 
found in both low-redshift and high-r edshift star-fo rming 
galaxies (e.g., [Taniguchi et al."2009l; lOesch et all 120101 : 
iGrazian et al.l l2012: Ono et al. 2012). For example, with 
a large sample pho to-selecte d LBG s at z 7 in the CAN- 
DELS fields, Gr azian et al.l poll ) found a well-defined 
relation, rhi oc L^/^. We illustrate the size-luminosity 
relation in our sample by displaying the best log-linear 
fit in Figure 2 (dashed lines). The best fitting results 
in the two panels are r\^\ oc L'^-^* and rhi. in ^ L'^'^, re- 
spectively. Although the re lation is obvious, the slopes 
are not as steep as that in IGrazian et al.l ()2012|) . The 
reason is that our galaxy sample is brighter. Our re- 
lation is derived from galaxies in the luminosity range 
of Afisoo < —19.5 mag, while the IGrazian et aP ()2012[ ) 
sample covers a range of Muv < — 18 mag. Their rela- 
tion largely depends on the galaxies fainter than -19.5 
mag, as seen in Figure 9 of their paper. In the brighter 
galaxies, rhi (or rhi, in) sho ws less of a trend with luminos- 
ity, as also pointed out bv IGrazian et al.l (|2012D . In fact, 
for the most luminous galaxies with M1500 —21 mag, 
their size shows little correlation with luminosity in our 
sample. These galaxies have likely reached a relatively 
mature structure even at z ~ 6. 

It should be pointed out that the size-luminosity re- 
lation is not significantly affected by our size measure- 
ments, i.e., we may lose more light at the edge of fainter 
galaxies, so that fainter galaxies tend to be smaller. This 
is because we only used high S/N detections (> Scr). 
Previous simulations and deep HST observations have 
shown tha t the above effect is minor for high S/N detec - 
tions (e.g.. lRavindranath et al]|2006l: lOesch et al.ll2010D . 
Furthermore, the size-luminosity relation in our sample is 



not caused by a selection effect, i.e., we may have missed 
some faint galaxies with large sizes. There are two rea- 
sons. One is that the galaxies were selected in optical 
images, they are bright in the optical broad bands (for 
LBGs) or narrow bands (for LAEs). The other one is 
that galaxies at the faint end tend to be very small, due 
to their hierarchical formation and evolution (Grazian et 
al. 2012; see also Windhorst et al. 2008 for detailed dis- 
cussion). Also note that the exclusion of 19 faint galaxies 
in our sample does not introduce bias to our results of 
the size-luminosity relation. Our HST near-IR data have 
relatively uniform depth (two orbits per band per point- 
ing), so the 8a cut indeed puts a flux limit on Afisoo in 
Figure 2, which does not affect our results in the bright 
region of Afisoo < —19.5 mag. 

3.2. Nonparametric Measurements of Morphology 

In this subsection, we will characterize galaxy structure 
and morphology using nonparametric methods, includ- 
ing the CAS (Concentration, Asymmetry, and Smooth- 
ness) system ([Consclicc 2003), the Gini coefficient G, and 
the M20 (second-order mom ent of the brightest 20% of 
a galaxy) parameters ( Lotz et al.ll2004[ ). These methods 
have been widely used for low-redshift galaxies. They 
usually provide reliable description of galaxy structure, 
and are able to distinguish different types of galaxies. To 
obtain accurate measurements of these quantities, two 
criteria are often required: high S/N ratios and large 
object sizes compared to the PSF size. For low-redshift 
galaxies, especially those in HST images, the sizes of 
galaxies are many times larger than PSF, so the two cri- 
teria are naturally met in deep HST images. At higher 
redshift, galaxies are fainter and smaller, so it is difficult 
to meet the two criteria, an d these parameters become 
less reliable. For example, IConselice fc Arnold! ()2009l ) 
tested the five parameters for galaxies at high redshift 
up to z ~ 6. They found that S fails to well describe 
dumpiness, while most of them can be used with cau- 
tion. 

3.2.1. CAS parameters 

Concentration (C) measures how compact the galaxy 
fight profile is. We adopt th e most recently and 
commonly use d definition (e.g., iBershadv et al.l [20001 : 
IConselicd[200l : 

C = 5loJ'-^), (1) 
V20J 

where rgo and r2o are the radii that contain 80% and 20% 
of the total galaxy fiux, respectively. Column 6 in Table 
1 shows the individual measurements of C. Compared 
to low-redshift galaxies, our galaxies have a relatively 
narrow range of C between ~2 and '^S. In particular, 
there is a lack of a population with C > 3, or a lack of 
highly-concentrated galaxies. This is likely because the 
measurements of rgo and r2o are not robust, due to the 
relatively low spatial resolution of the HST images for 
z > 6 galaxies. The inner radius r2o is usually smaller the 
PSF size. In particular, for highly-concentrated galaxies, 
r2o could be smaller than one pixel and is thus likely 
overestimated, so C is underestimated. Nevertheless, the 
derived C values generally refiect galaxy concentrations, 
as we will see in subsection 3.2.3. 
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Asymmet ry (A) measures how r otationally sym metric 
a galaxy is ()AbraIiam et al.lll996t IConselicel [20031 ) . It is 

calculated hy A = Agalaxy - ^background, Or: 



A 



El/ - I 



1801 



nB-B 



1801 



(2) 



In the first term ^galaxy , I is the galaxy image and /igo is 
the image rotated by 180° about the galaxy center. The 
galaxy center is determined to minimize ^galaxy The 
second term has the same meaning, except that it is for 
the sky background. A constant sky background has 
been subtracted from each galaxy image, so ^background 
is equal to zero in Equation 2. Usually, A is measured 
within a radius of 1.5 rp. For our galaxies, however, A 
is very sensitive to this radius due to the small sizes of 
the galaxies (or the small numbers of pixels occupied). 
We choose to measure A using the galaxy pixels assigned 
by segmentation maps given by SExtractor. Column 7 
in Table 1 gives the results. Because of our different 
method to calculate A, we cannot directly compare our 
results with others. We will show later that this A may 
still provide a reliable measurement of galaxy asymmetry. 

Smoothness {S) or dumpiness measures how clumpy a 
galaxy is (Conselice 2003). This parameter is sensitive to 
image resolution, and only works well for galaxies many 
times larger than the PSF size. Because our galaxies are 
either barely resolved or a few times the PSF size, we do 
not use S in this paper. 

3.2.2. Gmi and M20 parameters 

The Gini coefficient (G) measures how eve n the galaxy 
light distribution is ([Abraham et al.l [20031 : iLotz et all 
[2004D . It is calculated over the pixel s belonging to the 
galaxy. We use the formula given bv [Lotz et all (|200l 
to calculate G: 
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where is the absolute value of the i**^ pixel, and |/| 
is the average of the absolute values. Column 8 in Table 
1 shows the measurements. Compared to l ow redshift 
galax ies, e.g., z ~ 4 LBGs in the HUDF (jLotz et al.l 
[2006I) . our galaxies have slightly lower Gini coefficients. 
This may suggest that their light distrib ution is rela - 
tively even. Given the limitation of G fe.g. [Liskeii[2008[ ). 
an alternative explanation is that the image resolution 
is relatively poorer in our basically drizzled 2-orbit im- 
ages, compared to the well sampled and drizzled HUDF 
ACS images. So it cannot resolve subtle structures of the 
galaxies. 

The A/20 parameter, or the second-order moment of the 
brightest 20% of the galaxy, is similar to the concentra- 
tion C and measu res how the galaxy light is concentrated 
([Lotz et al.[ [20041 ). It is also calculated over the pixels 
belonging to the galaxy. The second-order moment of a 
pixel is defined by: 



Mi = Ii[(xi - XcY + iVi - VcT 



(4) 



where (xc, yc) is the galaxy center and (x^, yi) is the 
pixel position. We rank the pixels with decreasing order 
of flux, and select the brightest pixels whose total fiux is 
equal to 20% of the total galaxy flux. Then M20 is the 



total moment of the selected pixels normalized by the 
total moment of all the galaxy pixels: 

log (%^) , while ^ /, < 0.2/tot. (5) 
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The results are shown in Column 9 of Table 1 . Compared 
to low-redshift galaxies, our galaxies have a relatively 
narrow range of A/20, and is lack of highly-concentrated 
(very low A/20) galaxies. This is just like parameter G, 
and is caused by the low resolution of the images. We 
will see below the correlation between A/20 and G, as 
expected from their definitions. 

3.2.3. Relations among the morphological parameters 

As we mentioned earlier, the GAS system and the Gini 
and A/20 parameters are primarily used for low-redshift 
galaxies, and have not been used for z^Q galaxies. We 
for the first time used these morphological parameters 
for our z > Q galaxies. Most of our galaxies indeed do 
not meet the two criteria required to use these meth- 
ods, i.e., high S/N ratios and large galaxy sizes com- 
pared to the PSF size. Here we check the applicability 
of these parameters to z > 6 galaxies. It is meaningless 
to perform direct comparisons with low-redshift galax- 
ies, due to the much smaller galaxy sizes compared to 
the PSF in our images. The low resolution would sys- 
tematically bias the morphological measurements. As 
we have seen in the above subsections, G and A/20 were 
underestimated at least for hi ghly-concen t rated galax- 
ies. This was also noticed bv [Lotz et al.[ ([2006[ ). Our 
galaxies have systematically lower G values compared to 
low-redshift galaxies, likely because our images cannot 
resolve subtle structures. In addition, the definitions of 
A is slightly different for that in the literature. Therefore, 
instead of making comparison with low-redshift galaxies, 
we check these parameters internally, as most of the pa- 
rameters are expected to correlate with each other, and 
we have seen s uch correlations in low-re dshift galaxies 
(e.g. [Lotz et~ 2004, 2006; Cons ehce fc A rnold 2009). 

Figure 3 shows the relations among GA, G, A/20, and 
A/1500 for the galaxies in our sample. Smoothness S was 
not calculated because of the reason explained in subsec- 
tion 3.2.2. As expected, these parameters are correlated 
with each other: each parameter is correlated with one or 
more of the other parameters. There are several strong 
relations among them. The dashed lines in Figure 3 are 
the best linear fits to these relations, which are explained 
as follows. The moment A/20 is correlated well with G, A, 
and G. Both A/20 and G describe how the galaxy light is 
concentrated, and they are thus correlated by definition. 
The relations between A/20 and A k, G reflect that more 
concentrated galaxies have more rotationally symmetric 
profiles and more unevenly distributed light. The Gini 
coefficient G is correlated with G and A, in addition to 
A/20. Its relations with G and A indicate that galaxies 
with more concentrated or more rotationally symmetric 
profiles tend to have more unevenly distributed light. G 
is also correlated with A/1500. At higher luminosities, G 
has higher values. As we will see in the next section, 
many bright galaxies in our sample are interacting sys- 
tems, and thus show significantly uneven light distribu- 
tions. 

In Figure 3 we see all expected correlations among 
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these morphological parameters. Asymmetry A is de- 
fined in a slightly different way, but it is still correlated 
with M20 and G, indicating that our modified definition 
is reliable. Although it has been pointed out that con- 
centration C must been used with caution due to small 
galaxy sizes, C shows good relations with G and M20, 
and thus could be still reliable in most cases. We will 
see in the next subsection that the corrections applied to 
the above parameters are even better (smaller scatter) 
for bright galaxies in our sample. Therefore, morpholog- 
ical parameters such as C, A, G, and M20 could be still 
applicable for z > 6 galaxies, if used with caution. 

3.3. Interacting Systems 

One of the interesting morphological topics is to study 
interacting/merging systems, which traces hierarchical 
mass assembly in the CDM scenario. A close visual in- 
spection of Figure 1 shows that some galaxies are clearly 
extended with interacting or multi-component features. 
At lower redshift, these systems can be identified by 
pair coun ts fe.g.. iLe Fevre et al.l 120000 . the CAS sys- 
tem (e.g., 'Conselice"2003'), or the G and M20 parame- 
ters (e.g., Lotz ct al. 2008). Galaxies at 2; > 6 are faint 
and small, so it is difficult to properly distinguish reg- 
ular and interacting/merging systems. Although we de- 
rived morphological parameters GAGM20 in the above 
subsections, their ability to identify merging systems at 
z > 6 has never been examined. Therefore, we identify 
interacting/merging systems by visual inspection. Visual 
classification is the earliest way for galaxy classification, 
and in many cases is still the best way to identify merging 
systems at high redshift. 

We consider the following two types of galaxies as 
candidate interacting/merging systems: 1) galaxies with 
two or more distinct cores; and 2) galaxies with ex- 
tended/elongated features and/or long tails. We iden- 
tify these systems in 25 relatively bright (M1500 < —20.5 
mag) galaxies. For fainter galaxies, our images are not 
deep enough to properly identify all faint components 
or extended features due to low S/N ratios. Figures 4 
and 5 show 12 galaxies that are identified as interacting 
systems with M1500 < —20.5 mag. Figure 4 shows 11 
galaxies in the stacked images, and Figure 5 shows the 
z = 6.96 LAE (no. 62) in the two J bands. The red pro- 
files are the SB contours of the rest-frame UV emission. 
Six of them clear show double or multiple clumps, in- 
cluding no. 4, 24, 34, 49, 62, and 67. They usually have 
one bright core and one or more fainter clumps. No. 62 
and 67 are particularly interesting. No. 62 in Figure 5 
has almost two identical components, and no. 67 has 
three widely-separated cores. More discussion is given in 
the next subsection. Another 3 galaxies, including no. 
36, 58, and 61, do not clearly show multiple clumps, but 
have long tails like tidal tails seen in low-redshift merg- 
ing galaxies. The rest of the 12 systems (no. 15, 44, and 
47) do not show multiple components or tails, but they 
are rather extended and elongated. They could be in the 
end of the merging process. 

We estimate the fraction of mergers among galaxies 
with M1500 < -20.5 mag. The fraction is 48%, or 36% 
if we exclude the three galaxies that do not show mul- 
tiple components or tails. The fraction is even higher 
in the galaxies with M1500 < —21 mag. We have 18 
galaxies in this magnitude range, and 10 of them are 



mergers. The fraction of mergers is 56%, or 39% if we 
exclude the three galaxies mentioned above. This is con- 
sistent with the fractions in the brightest galaxies at low 
redshift of z ~ 2 — 3. For example, the fractio n of merg- 
ers in the Mb < —21 mag galaxy sample of Conselic^ 
(|200l is 40%-50%. The typical fraction i n the brightest 
galaxies in the sample of lLaw et al.l ()2012[ ) is also ^ 40%. 
The merger fraction in fainter galaxies is smaller, because 
merger systems consist of multiple components, and have 
stronger SFRs and UV emission. For the galaxies fainter 
than M1500 — —20.5 mag in our sample, our images are 
deep enough to identify double-core systems with com- 
parable emission. We find that these systems are very 
rare at M1500 — —20.5 mag. Studies in deeper fields 
clearly indicate a low mer ger fraction in low-luminosity 
galaxies. For example, Oes ch et al.l (j2010D presented the 
morphology of 16 z > 7 LBGs in the HUDF, and only 
found two galaxies with extended features. 

3.3.1. Notes on individual objects 

Galaxy no. 62. No. 62 is a z = 6.964 LAE. It i s 
the first spec-confirmed LAE at z ~ 7 (|Ive et al.ll2006f) . 
Figure 5 shows that it has two similar components in 
the both J bands. We use the two -dimensional fitting 
algorithm GALFIT (|Peng et al.ll2002f ) to model the two 
components. Two Sersic functions are fitted to the two 
components (two left-hand images in Figure 5) simulta- 
neously. The middle images in Figure 5 show the best 
model fits, and the residuals are on the right-hand side. 
The two components in the both bands can be well de- 
scribed by the Sersic function. The separatio n between 
the t wo cores is about 0'.'2 (^ 1 kpc) (see also lCai et al.l 

[mil). 

Galaxy no. 61. No. 67 is a LAE at z = 6 .595. It was 
discovered as a giant LAE bv lOuchi et al.l (|2009[ ). It is 
one of the brightest galaxies in terms of both Lya lumi- 
nosity and UV continuum luminosity. The most striking 
feature is the three well separated cores lined up. The 
central core is relatively weak. The brighter side cores 
are even resolved in our ground-based z— and y-band 
images. The separate between the two side cores reaches 
~ 1'.'2, or 7 kpc. This is the largest separation we have 
seen at 2 > 6. This object also has strong emission in the 
IRAC band s, suggesting a large stellar mass of > IO^^Mq 
(jOuchi et a l. 2009). 

3.3.2. Morpfiology of interacting systems 

As we mentioned, the above interacting galaxies were 
identified with visual inspection. Morphological param- 
eters were not used because the criteria of interacting 
systems at z > 6 are unclear. We check our interacting 
systems in the parameter space in Figure 6. Figure 6 
is similar to Figure 4, but only plots the galaxies with 
-^1500 < —20.5 mag in our sample. The red squares rep- 
resent the interacting galaxies. It is clear in Figure 6 that 
the interacting galaxies are almost indistinguishable from 
the rest in the space of all parameters except A. Asym- 
metry A is the only effective parameter to identify inter- 
acting systems. The galaxies with the largest A in our 
sample are nearly all interacting systems. A is widely 
accepted as the mo st efficient way to identi fy mergers 
at low redshift (e.g. lConselicell2003l:IC onsclice fc Arnoldl 
l2009HLotz eralll2006ll2008HLaw et al.li2012)) . The typi- 
cal criterion for low-redshift mergers is A > 0.30 — 0.35. 
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Although all the mergers in our sample have A values 
greater than 0.35, the criterion also includes a large num- 
ber of regular galaxies. This is not a good criterion for 
our galaxy sample, because of the slightly different defi- 
nition of A and the low spatial resolution of our images 
as we emphasized earlier. A better criterion in Figure 6 
is probably A > 0.45, which includes the majority of the 
mergers, but few of non-merging galaxies. 

Another striking feature in Figure 6 is the much better 
correlations among the morphological parameters com- 
pared to Figure 4. The scatter in the relations is much 



smaller for these bright galaxies with M; 



1500 



< 



-20.5 



mag. There are two explanations. One is that bright 
galaxies have higher S /N ratios and cover more pixels in 
the HST images, so their morphological measurements 
are more robust. The other one is that galaxies with 
different luminosities occupy slightly different parameter 
space, so a sample covering a smaller luminosity range 
shows a smaller scatter in Figure 6. The real reason 
is very likely the combination of the two. In any case. 
Figure 6 strengthens our earlier conclusion that morpho- 
logical parameters CAGM20 could be still applicable for 
z > 6 galaxies. 

4. Lja MORPHOLOGY 

In this section we will study the Lya morphology of 
LAEs using our ground-based narrow-band images. We 
will not measure structure and morphology for individ- 
ual galaxies, since they are mostly point-like objects in 
the ground-based images. Although these images have 
excellent PSF sizes of ~ 0'.'5 - 0'.'7, the PSF sizes are stiU 
much larger those of the HST images. Therefore, we 
will focus on Lya halos around LAEs, which could ex- 
tend many arcsec from the objects. We will also compare 
the positions of Lya emission with those of UV contin- 
uum emission, and find any possible positional difference 
between the two. 

4.1. Lya Halos 

Because of the resonant scattering of Lya photons by 
neutral hydrogen, Lya emission could form large diffuse 
Lya h alos around high-redshift galaxies. S teidel et al.l 
(poll first found very extended Lya halos in a sample 
of luminous galaxies at 2 < z < 3. The galaxies were UV 
continuum-selected, but more than a half of them show 
net Lya emission and ~ 20% have Lya EW greater than 
20 A. They were able to find large Lya emission halos 
(> 80 kpc) in the stacked images of all sub-samples of 
their galaxies. They further claimed that all LBGs would 
be classified as LAEs or Lya blo bs, if imaging data are 
deep enough to detect Lya halos. iMatsuda et al.l (|2012[ ) 
confirmed the existence of extended Lya halos around 
z ~ 3 galaxies. They used more than 2000 LAEs at 
z ~ 3.1, and grouped them into sub-samples based on lu- 
minosity and surface overdensity. They stacked narrow- 
band (Lya) images for each sub-sample, and found that 
all stacked images sh ow extended (> 60 kpc ) Lya emis- 
sion halos. Recently iFeldmeier et al.l ()2013D found that 
the existence of Lya halos around LAEs is not convinc- 
ing. They also used a large sample of a few hundred 
LAEs at z ~ 2.1 and 3.1. They payed particular at- 
tention to systematic effects from large-radius PSF and 
large-scale flat fielding, etc. When these effects were 
taken into account, they did not find strong evidence of 



extended Lya halos in the stacked narrow-band images 
at either redshift. They tried a few ways to reconcile the 
discrepancy between their results and the previous re- 
sults, yet the reason of the discrepancy is still not clear. 

Stacking of narrow-band images has not been done for 
z > 6 galaxies. On the other hand, cosmological sim- 
ulations have predicted the existenc e of extended Lya 
emission around z ^ 6 galaxies (e.g.JZheng et al. 120111 : 
iDiikstra fc Krameiil2012tlJeeson-Daniel et al.ll2012D . For 



example, by including the resonant scattering of Lya 
photons in bot h circumgalactic me dia and intergalactic 
medium fIGM). IZheng et al.l(|2011h showed that the Lya 
emitting halo in a high-redshift galaxy can extend up to 
1 Mpc. They further pointed out that such halos could 
be detected by stacking 100 z ~ 5.7 LAEs in 4-hr expo- 
sure narrow-band images in the SXDS. Here we combine 
the narrow-band (Lya) images of LAEs at z c± 5.7 and 
6.5. 

In order to detect diffuse Lya halos around LAEs, 
we make use of all the known (spec-confirmed) LAEs 
at z ~ 5.7 and 6.5 in th e SD F and SXDS fields 
from iKashikawa et all (|2011l ) and lOuchi et al.l (|2008l 
l20Tnh . The narrow-band images were taken with the 
Subaru telescope. The total integration time in the 
NB81 6 and NB921 bands a re 10 and 15 hr for the 
SDF (IKashikawa et al.ll2004f) . and -4 hr and -10 hr 
depth slig htly varies among five pointings) for the SXDS 
Ouchi et al. 2008, 2010). The data reduction of the SDF 
images were presented in Section 2 of Paper I. The re- 
duction of the SXDS images were done in the same way. 
These narrow-band images have great depth with excel- 
lent PSF sizes of 0'.'5 ~ 0'.'7. We reject a small number of 
galaxies that are either very faint or blended with nearby 
bright objects. The final sample contains 43 LAEs at 
z ~ 5.7 and 40 LAEs at z ~ 6.5. 

To stack these images, we first cut image stamps for all 
individual objects. We then re-sample image stamps, so 
that the objects are in the centers of each image. After 
all other objects in the images are masked out, we co-add 
(average) the images with sigma-clipping {5a rejection). 
Figure 7 shows the stacked narrow-band Lya images at 
the two redshifts. We did not subtract away continuum 
images. As we showed in Paper I, our LAEs have Lya 
EWs greater than 20 A. The median EW value is 74 A, so 
their narrow-band photometry is completely dominated 
by the Lya emission. Subtracting continuum images 
would significantly increase the background noise. In the 
above procedure we did not scale the objects to the same 
magnitude either, since we do not know how the SB of 
the Lya halos scales with the total Lya flux in a LAE. We 
perform a series of tests by stacking the images in differ- 
ent ways, including image stacking with median, with ob- 
jects scaled to the same magnitude, and with bright ob- 
jects only. The results are all very similar to the images 
shown in Figure 7. Our stacked images of LAEs reach a 
depth of roughly 1.2 x 10~^^ erg s~^ cm~^ arcsec"^ {la) 
in both bands (see also Figure 8). We should point out 
that this depth is shallower than what we predicted as- 
suming that the noise goes down with the square root 
of the image number. The reason is that the input im- 
ages are not purely blank images, whose noise is dom- 
inated by background Poisson noise. These images are 
extremely deep and crowded, so LAEs have many nearby 
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objects. Although we were able to mask out the bright 
pixels of the nearby objects, we could not remove the 
light outside these pixels that may extend many pixels 
away. When we combine many images, these nearby ob- 
jects contribute significant unavoidable no ise to the co- 
added images (e.g.. Windhors t et al.ll2008l ). 

Figure 7 clearly shows that the stacked Lya emission 
is compact at the both redshifts, and do not show very 
extended diffuse halos. As a reference, we also show in 
Figure 7 the stacked images of stars (point sources) . The 
stars are chosen to be separated, bright, but not satu- 
rated, in the same narrow-band images. They are com- 
bined in the same way as we did for LAEs. Figure 8 
shows the radial profiles of the stacked LAEs (solid lines 
with error bars) as well as the profiles of the stars (dotted 
lines). The Lya profiles roughly follow the PSF profiles, 
and their SB reach zero at r > 2". In the upper panel, 
the profile FWHM of the stacked z ~ 5.7 LAE and star 
are 0'.'49 and 0'.'62, respectively. In the lower panel the 
two FWHM values are 0'.'61 for the z ~ 6.5 LAE and 
0'.'77 for the star. The LAE profiles are broader than 
the PSF sizes by ^ 26%, and exhibit slightly longer tails 
than the PSF profiles do. This is simply because galaxies 
are not point sources, and indicates that the Lya emis- 
sion is resolved but not very extended. In the previous 
sections we show that our L AEs have a range o f sizes in 
their rest-frame UV emission. iHathi et all (|2008t) stacked 
broad-band HUDF images for z 4 — 6 LBGs and found 
that the galaxy SB profiles are apparently broader than 
the PSF profiles. Finkelstein et al.l (pOll.) found that 
two z ~ 4.4 LAEs have large sizes in Lya than in UV 
continuum. So the slightly broader Lya SB profiles com- 
pared to the PSF are just the nature of galaxies, and are 
not likely caused by diffuse Lya halos predicted or found 
in previous studies. The z — 5.7 profile seems to exhibit 
slightly more extended radius than its PSF profile. Given 
the 1(7 limit of ^ 1.2 x 10~^^ ergs~^ cm~^ arcsec"^, its 
SB is still consistent with zero at r > 2". Therefore, 
Figure 8 does not show convincing evidence of extended 
Lya halos. 

It is difficult to answer whether our stacked images are 
deep enough to detect Lya halos if they do exist at the 
two redshifts. Our images are certainly de ep enough to 
detect the LBG halos at z ~ 3 reported by 'Steide l et al. 
(120111) and the LAE halos at z ~ 3.1 in 'Matsud a et al. 
l20m . The depth to detect z > 6 LAE halos is ob- 
servationally unkn own. For cosmological simulations, 
iZheng et al.l (|20lH ) predicted Lya halo sizes in z = 5.7 
LAEs. They found two characteristic scales for the halos. 
The inner steeper one extends to 3" — 4", and the outer 
flatter one extends to a few tens of arcsec. While our im- 
ages are not deep enough to detect the outer halos, they 
are almost deep eno ugh to detect the in ner scale halos as 
seen in Figure 3 of IZheng et al.l (|2011| ). where the Lya 
radial profiles are shown for LA Es in dark matter halos 
of ~ IO^Mq. As pointed out bv lZheng et"all (|20lH) . the 
size of diffuse Lya emission also depends on the mass of 
dark matter halo. If the average mass of the dark matter 
halos in our LAEs is smaller than 10^^ Mq, our current 
data may not be able to detect the diffuse Lya emission. 

It is also likely that the Lya halos (if exist) have been 
diluted to a much lower level during the construction 
of the stacked images. The stacked images can prop- 
erly recover Lya halos only when halos are smoothly and 



symmetrically distributed around galaxies. If halos are 
highly asymmetric and/or clumpy, the emission of halos 
will be significantly diluted in average stacked images, 
and could totally disappear in median stacked images. 
From the observations of z ~ 6 quasars or cosmological 
simulations, we know that t he distribution of IGM at 
z ~ 6 is inhomogeneous (e.g. iFan et al.l 120061: iMesingeii 
|2010[) . If the distribution of IGM affects the shape of 
Lya halos (via resonant scattering), the distribution of 
Lya halos is also likely asymmetric and clumpy. 

Finally, it is possible that these LAEs do not have 
extended Lya halos, or t hat their ha l o emi ssion is not 
as strong as predicted bv lZheng et al.l (1201 ID. especially 
when dust is taken into account (e.g. iFinkelstein et al.l 
120 111) . IZheng et all ()2011l) did not consider dust in their 
simulations. We know that high-redshift LAEs are not 
free of dust. In particular, the brightest galaxies may ex- 
hibit significant dust extinction, as implied by their UV 
colors (see Paper I). When Lya photons are resonantly 
scattered by dusty neu tral hydrogen, the L ya emission is 
substantially reduced ()Yaiima et al.ll20l"2l ). The reduc- 
tion is more severe at larger distance from the object, 
because photons at larger distance need to pass through 
more dust before they escape. This process would signif- 
icantly reduce the visibility of possible diffuse Lya emis- 
sion, and makes it much more difficult to detect it. A 
much larger sample of LAEs is needed to answer this 
question. 

4.2. Lya- Continuum Misalignment 

The comparison between the positions of UV contin- 
uum and Lya emission provides useful information on 
how Lya photons escape from a galaxy. Lya and UV 
continuum photons are usually come from the same star- 
forming regions, although Lya photons are likely more 
sensitive to the regions with more recent star-forming 
activity. As we mentioned earlier, Lya emission is com- 
plicated by resonant scattering and IGM absorption. So 
the observed position of Lya emission could be different 
from the position of UV continuum emission. For ex- 
ample, a large positional difference has been found in a 
z = 3.334 galaxy (?). Due to the small sizes of high- 
redshift galaxies, current ground-based observations are 
not able to detect these positional differences. We rely 
on HST, which has observed a large number of high- 
redshift galaxies. HST observations were mostly made 
for rest-frame UV continuum emission, and there is usu- 
ally no suitable HST narrow-band filters for Lya emis- 
sion. One ex ample for the HST imag ing of Lya emission 
is the work of lFinkelstein et all ()2011l ) , who observed the 
Lya emission of a small sample of ~ 4.4 LAEs with a 
narrow-band filter. They did not find strong evidence of 
positional misalignment between UV and Lya emission. 

We use our large sample of LAEs to search for possi- 
ble positional offsets between UV and Lya emission at 
z ~ 5.7 and 6.5. We use our HST images as UV contin- 
uum images and Subaru narrow-band images as Lya im- 
ages. As mentioned above, the narrow-band images have 
excellent PSF FWHM sizes around 0'.'5 - 0'.'7. During 
the construction of the HST images (Paper I), we have 
matched the coordinates of the HST images to those of 
the optical images. To avoid any large-scale variation, 
we refine the coordinates of the HST images. For each 
LAE, we find 10-20 nearby objects that are relatively 
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bright and round. We then match the positions of the 
nearby objects in the two sets of images. The typical 
refinement is smaller than the size of one pixel (0'.'06). 
The uncertainty in the object positions, derived from the 
distribution of the nearby objects, is about the size of 1- 
2 pixels. We plot the Lya positions on top of the UV 
continuum positions. Figure 9 shows a few examples of 
bright galaxies. The red profiles are the contours of the 
UV emission SB, and the green crosses indicate the po- 
sitions (and la uncertainties) of the Lya emission. 

For the majority of the LAEs in our sample, the posi- 
tions of the UV continuum and Lya emission agree with 
each other. In particular, we do not find positional mis- 
alignment at a significance level of > 2a among almost 
all compact and round LAEs. Object no. 3 in Figure 
9 is a typical example, in which the center of the Lya 
emission is close to the center of the continuum emis- 
sion. For the merging/interacting systems, however, we 
see significant positional differences. Figure 9 shows the 
examples of these systems. They exhibit a variety of Lya 
positions relative to the peak positionfH of the UV con- 
tinuum emission: 1) Lya positions close to the positions 
of the brightest components in the UV images, including 
galaxies no. 15, 47, and 61; 2) Lya positions close to the 
positions of the fainter components or merger tails in the 
UV images, including no. 14, 44, and 58; 3) Lya posi- 
tions are somewhere between the positions of the bright 
and faint components, but closer to the bright compo- 
nents, including no. 4 and 49. No. 62 and 67 are again 
interesting. The Lya position of no. 62 is roughly in 
the middle of the two similarly bright components. No. 
67 does not show three distinct Lya emission clumps, as 
its UV emission does. Instead, it shows a single bright 
Lya emission core with some extended features. The 
Lya center is not at any of the three UV clumps. It is 
between the left and central clumps, and slightly closer 
to the central one. 

Our results suggest that in compact galaxies the ob- 
served location of Lya emission does not deviate from 
its original position, while in merging/interacting sys- 
tems, the observed Lya location could be significantly 
different from its original position, without preferential 
positions of offsets. If the final location of Lya emis- 
sion is determined by the process of resonant scattering, 
our results can be explai ned by the interstellar medium 
(ISM) distribution (e.g. iFink elstein et all 1201 1( ). In a 
non-disturbed galaxy, the ISM distribution is relatively 
symmetric around the object (not necessarily uniform, 
it could be clumpy). The random scattering of photons 
does not have preferred directions, so the observed loca- 
tion is still close to its original location. In an interacting 
system, the ISM is re-distributed by merging activity. 
The distribution of the disturbed ISM is therefore not 
symmetric any more. This results in an offset of the 
observed Lya position. 

This result could have an impac t on spectr o scopy 
of bright LBGs at z > 6 (e.g. iStark et al.l 120111: 
iCurtis-Lake et all 120121 ). While most of the positional 

For an interacting galaxy, the position of its peak emission 
(as seen from the SB contour) could be very different from the 
position of the overall galaxy emission measured by, for example, 
fitting Gaussian to the marginal x,y distributions (used by IRAF 
DAOPHOT). Here our positions refer to the positions of peak emis- 
sion. 



offsets in our sample are smaller than 0'.'2, at least two 
are around 0'.'3 — 0'.'4. If one uses a < 1" slit to identify 
their Lya emission lines, based on the positions of con- 
tinuum emission, one could miss them due to the large 
offsets. However, glaxies with such large offsets are very 
rare, so this result will not largely reduce the success rate 
of identifying Lya lines in LBGs. 

5. DISCUSSION AND SUMMARY 

The comparison in Section 3.1 shows that the sizes of 
our galaxies are consistent with those of bright z > 6 
LAEs and LBGs in previous studies. Th is disagrees with 
the claim of Dow-Hv gelun d et al.l (|2007[ ) that LAEs are 
more compact than LBGs. This is likely because the 
iDow-Hvgelund et al.l ()2007| ) sample is very small while 
LAEs have a large range of sizes. It is indeed difficult to 
make proper comparisons without a large sample, as the 
galaxy size depends on redshift and luminosity. Most 
of our galaxies are LAEs. An other large LAE sampl e 
is the z ~ 5.7 LAE sample by i Taniguchi et al.l (|2009| 1. 
who observed a number of LAEs at z ~ 5.7 with the 
HST ACS F814W fiffer. Ahhough this filter includes 
Lya emission, it is so wide that its emission is dominated 
by UV continuum. They found that the average intrinsic 
size of the sample is 0'.'13, the same as ours. Among 
many photo-selected LBG s amples at z ^. 6, a recent 
large sample is the sample of lGrazian et al.l ()2012| ). This 
sample contains a number of bright (as well as faint), 
photo-selected LBGs at 2; ~ 6 - 7 in the CANDELS field. 
The sizes of these galaxies are well consistent with those 
of our galaxies at the same luminosities. Therefore, we 
conclude that LAEs and LBGs with similar luminosities 
have similar physical sizes. 

In this paper (and Paper I), galaxies found by the 
narrow-band technique are defined as LAEs and those 
found by the dropout technique are defined as LBGs. As 
we discussed in Paper I, this classification only refiects 
the methodology that we apply to select galaxies. It does 
not mean that the two types of galaxies are intrinsically 
different. In Section 5.3 of Paper I, when we derived the 
UV continuum luminosity function of LAEs, we used an- 
other popular definition of LAEs based on the Lya EW, 
i.e., a galaxy is a LAE if its Lya EW is greater than 20 A. 
With this definition, almost all the galaxies in our sample 
are LAEs. This definition is physically more meaningful, 
but observationally difficult, because one can easily ob- 
tain a flux-limited sample, not a EW-limited sample. We 
have 16 LBGs (former definition) in our sample. They 
are not typical LBGs. They are spec-confirmed, and thus 
only represent those with strong Lya emission. In Paper 
I we found that our LAEs and LBGs are indistinguishable 
in many aspects of the Lya and UV continuum proper- 
ties. In this paper, we further found that these LAEs 
and LBGs have similar physical sizes and morphological 
parameters. This confirms one of our conclusions in Pa- 
per I that LAEs are a subset of LBGs with strong Lya 
emission lines. 

This paper is the second in a series presenting the 
physical properties of a large sample of spec-confirmed 
galaxies at z > 6. The sample consists of 51 LAEs and 
16 LBGs, and represents the most luminous galaxies in 
terms of Lya luminosity (for LAEs) or UV continuum 
luminosity (for LBGs) in this redshift range. In Paper 
I we derived various properties of rest-frame UV con- 
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tinuum and Lya emission. In this paper we have con- 
ducted a detailed structural and morphological study of 
the galaxies using deep HST near-IR images and Sub- 
aru narrow-band images. In order to measure the mor- 
phology of rest-frame UV continuum emission, we con- 
structed a stacked HST image for each object by combin- 
ing its J- and i?-band images. UV morphology was then 
measured for those with > 8ct detections in the stacked 
HST images. We used half-light radius rhi to describe 
the sizes of galaxies. The observed sizes of our galaxies, 
without correction for PSF broadening, are between O'.'l 
and O'.'S (or 0.6-1.6 kpc), with a median value of 0'.'16. 
The intrinsic sizes rhi^n in our sample are from < 0'.'05 
(< 0.3 kpc) to ^ 0'.'26 1.5 kpc), with a median value 
of O'.'IS (~0.75 kpc). These values are consistent with 
those of bright, photo-selected LBGs at similar rcdshifts. 
Additionally, more luminous galaxies tend to have larger 
sizes, exhibiting a size-luminosity relation rhi^n k L^'^. 
The slope 0.2 is significantly flatter than those in fainter 
LBG samples. Our objects show a wide range of mor- 
phology in the HST images, including compact galaxies 
and double/multiple component systems. The bright- 
est galaxies in the sample have a large fraction of merg- 
ing/interacting systems. The fraction of mergers reaches 
40 — 50% at Mi5oo < —20.5 mag, suggesting hierarchical 
build-up of galaxies at high redshift. 

We for the first time described the structure and mor- 
phology of z > 6 galaxies using nonparametric methods, 
including the CAS system, the Gini and M20 parame- 
ters. Compared to low-redshift galaxies, z > 6 galaxies 
appear faint and small in HST images, so they usually 
do not meet the criteria required by these nonparamet- 
ric methods for low-redshift galaxies. Because of this 
reason, our galaxies are located in slightly different re- 
gions in the parameter space. However, we found strong 
correlations between these parameters, as expected from 
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their definitions. This suggests the nonparametric meth- 
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only effective parameter to identify merging systems. 

We for the first time searched for Lya emission halos 
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bined a large number of narrow-band images for LAEs 
at z ~ 5.7 and 6.5, respectively. The stacked images 
reached a depth of ^ 1.2 x 10""'^^ ergs~^ cm~^ arcsec"^ 
(Ict). We did not find evidence of extended diffuse Lya 
emission as predicted by cosmological simulations. It is 
possible that our images are still not deep enough to 
detect Lya emission halos, or that the Lya halo emis- 
sion has been diluted to a much lower level during the 
construction of the stacked images. It is also possible 
that the halo emission is not as strong as predicted. A 
much larger LAE sample is needed to solve this ques- 
tion. We also investigated positional differences between 
the rest-frame UV continuum emission and Lya emission 
in LAEs, using the HST images and optical narrow-band 
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of offsets. It was explained by the distribution of the 
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Morphological Parameters of the Galaxies 
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+27:24:52.04 


6.556 


0.196 


2.25 


0.44 


0.35 


-1.18 


58 


13:24:55.772 


+27:40:15.31 


6.583 


0.233 


2.57 


0.54 


0.44 


-1.18 


59 


13:24:58.508 


+27:39:12.92 


6.585 


0.208 


2.86 


0.46 


0.40 


-1.03 


61 


13:23:49.186 


+27:32:11.39 


6.599 


0.183 


2.10 


0.49 


0.38 


-1.37 


62 


13:23:53.054 


+27:16:30.75 


6.964 


0.179 


2.25 


0.29 


0.57 


-1.15 


63 


13:24:15.678 


+27:30:57.79 


6.023 


0.144 


2.36 


0.25 


0.50 


-1.34 


64 


13:24:40.239 


+27:25:53.11 


6.116 


0.141 


2.44 


0.36 


0.51 


-1.57 


66 


13:23:52.680 


+27:16:21.76 


6.575 


0.173 


2.34 


0.45 


0.34 


-1.04 


67 


13:24:10.817 


+27:19:28.08 


6.595 


0.155 


2.17 


0.60 


0.40 


-0.69 



Note. — The half-light radii rhi are in units of arcsec. 
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Fig. 1. — Thumbnail images of the 47 (out of 67) galaxies that have more than 8cr detections in the stacked HST images, shown in 
order of increasing redshift. The image size is 25 pixel square, which corresponds to 2'.'5 X 2'.'5 in NICMOS images (no. 25, 30, and 31), or 
1'.'5 X 1^5 in WFC3 images (the rest). 
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Fig. 2. — Half-light radius as a function UV continuum luminosity Mi^qq. The upper panel shows the observed radius rj^i without 
correction for PSF broadening. The dotted line indicates the PSF size in our HST WFC3 images. The lower panel shows the intrinsic 
radius rhi,in after correction for PSF broadening. The blue and red circles represent the LAEs at 2 ~ 5.7 and 6.5 (including z ~ 7), 
respectively, and the green circles represent the LBGs at 2 ~ 6. The star represents the stacked object of 19 faint galaxies. The dashed lines 
are the best log-linear fits to all the data points, suggesting a good size-luminosity relation: brighter objects tend to have larger physical 
sizes. 
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Fig. 3. — Correlations among C, A, G, M20, and M1500 for the galaxies in our sample. Each parameter is correlated with one or more 
of the other parameters. There are also several strong relations among them. The dashed lines are the best linear fits to the relations. We 
see all expected correlations among these parameters, suggesting that these parameters could be still applicable for z > 6 galaxies, if used 
with caution. 
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Fig. 4. — Interacting/merging galaxies in the HST images. The image size is 1'.'26 X 1'.'26. The red profiles are the contours of SB. Each 
contour starts at 85% of the peak value with an interval of 20%. These galaxies have extended and elongated features, and/or have double 
or multiple components. 
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Fig. 5. — Images of the z = 6.96 LAE (no. 62). The two left-most images show the galaxy in the two J bands (Jiio and J125). The two 
images in the middle are the best-fit model (two Sersic functions for the two components) galaxies. The residuals are on the right-hand 
side. The separation between the two components is about 0''2 (~ 1 kpc). 
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Fig. 6. — Same as Figure 4, but for galaxies with A/1500 < —20.5 mag. The red squares represent the interacting galaxies identified in 
our sample. The figure also shows better corrections (smaller scatter) among the parameters compared to Figure 4. 



18 




Fig. 7. — Stacked narrow-band images for LAEs and stars (or point sources) in two bands NB816 and NB921. The image size is 20" X 20". 
The images have the same intensity scale. The PSF FWHM sizes derived from the two stacked stars arc 0"49 and 0"61, respectively. The 
FWHM of the stacked LAEs are 0'.'61 and 0'.'77, respectively. They are larger than the PSF sizes by 26%, because LAEs are not point 
sources. The stacked LAEs do not show diflfuse Lya halos. 
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Fig. 8. — Radial profiles of the stacked images in Figure 7. The units of SB is lO^-''^ ergs^^ cm^"^ arcsec^"^. The solid profiles with Ic 
error bars represent the stacked LAEs, and the dashed profiles represent the stacked stars. The insets show the radial profiles on a log scale. 
The LAE profiles are broader than the PSF sizes, and exhibit slightly longer tails than the PSF profiles do, meaning that the Lya emission 
is resolved. This is because galaxies are not point sources, as explained in Section 4.1. At r > 2", the SB in both cases is consistent with 
zero within Icr errors. We do not see very extended halos of Lya emission. 
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Fig. 9. — Positions of UV continuum and hya emission for a sample of bright galaxies. The red contours and their scales are similar to 
those shown in Figure 4. They display the UV continuum emission seen in the HST images. The green crosses indicate the positions (and 
Icr uncertainties) of the Lyo emission from our ground-based narrow-band images. The first object no. 3 represents a typical compact 
galaxy, whose positions of UV continuum and Lya emission agree with each other. The rest of the objects are merging/interacting systems, 
which show a variety of hya position offsets relative to the UV continuum positions, including significant positional misalignment. 



